AD-A216  966 


NIC  FILE  COM 


Technical  Document  1687 
November  1989 

Development  of  a 
Progressive  Failure  Model 
for  Strength  of  Laminated 
Composite  Structure 


P.  Y.  Tang 


Approved  for  public 


DTIC 

ELECTE 
JAN  2  4 1990 


release;  distribution  Is  unlimited. 


90  01  2S  IS 


NAVAL  OCEAN  SYSTEMS  CENTER 

San  Diego,  California  92152-5000 

J.  D.  FONTANA,  CAPT,  USN  R.  M.  HILLYER 

Commander  Technical  Director 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  by  the  Structural  Mechanical  Analysis  and  Design  Branch, 
Code  931,  Naval  Ocean  Systems  Center,  under  Independent  Exploratory  Development  Pro¬ 
grams,  OCNR-20T,  Arlington,  VA 


Released  by 
P.  A.  Jung,  Head 
Structural  Mechanical  Analysis 
and  Design  Branch 


Under  authority  of 
C.  L.  Ward,  Jr.,  Head 
Design  and  Development 
Division 


JG 


REPORT  DOCUMENTATION  PAGE 


form  tpprwmi 
OMB  No,  0704-01 W 


-  3I^  opening  boidu  foe  tn*  oculoalon  d  Htam«flon  k  mttmstn  to  muqi  1  how  pw  iwpow.  InctoUng  o»  tlmo  tof  nZySHHnq  mwucBon*.  tUSSStg  mSStg  «Sa  t 


I.  AGENCY  USE  ONLY  |llM  bl*n*J 


4.  TITLE  AND  SUBTITLE 


04-01881.  W 


2.  REPORT  DATE 

November  1989 


jj  fcjp 


3.  REPOm  TYPE  ANO  OATES  COWERED 

Interim 


5.  FUNONQ  NUMUERS 


DEVELOPMENT  OF  A  PROGRESSIVE  FAILURE  MODEL  FOR  STRENGTH 
OF  LAMINATED  COMPOSITE  STRUCTURE 


8.  ALT'  Ofl(S) 

D.  Y.  Tang 


7.  r  ERFORMING  ORGANIZATION  NAME(S)  ANO  ADORESS<ES) 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152-5000 


r  XJNSORING/MONITORING  AGENCY  NAME(S)  ANO  ADORESS<ES) 

Independent  Exploratory  Development  Programs,  OCNR-20T 
Arlington,  VA  22217 


PE:  602936N 
WU:  RV36I21 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NOSCTD  1687 


10.  SPONSORINO/MONfTORlNO 
AGENCY  REPORT  NUMBER 


1 2a.  01. .  fRlSUmON/ AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Uvomcrn  200  worts) 

^~D^This  report  briefly  describes  the  background  information  and  current  work  plan  for  a  multiyear  effort  to  develop  a 
j  ogressive  failure  model  for  predicting  strength  of  laminated  composite  structure  under  static  loads  and  fractured  in  the 
i  -plane  failure  modes.  It  also  highlights  the  intermediate  results  obtained  in  FY89.  , ..  iir,r.rl<r  i  <rr  \  . 

4 


14.  SUBJECT  TERMS 

laminated  composite  structural  strength 
progressive  failure  model 

piecewise  quadratic  strength  tensor  failure  criterion 
stiffness  degradation  model _ 


1 7.  SECURITY  CLASSIFICATION  1 8.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE 

UNCLASSIFIED  UNCLASSIFIED 


IS.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


IS  NUMBER  OF  PAGES 


19.  PWCE  COTE 


M  UAKTAOON  OF  ABSTRACT 


Unlimited 


NON  7540-01-280-5500 


S5d.\Vkvi\'in'  M 


CONTENTS 


1.0  INTRODUCTION .  1 

2.0  BACKGROUND  .  2 

2.1  Composite  Stress  Analysis .  2 

2.2  The  Piecewise  Quadratic  Strength  Tensor  Failure  Criterion .  4 

2.3  Progressive  Failure  Model .  5 

3.0  DEVELOPMENT  PLAN .  7 

3.1  Work  Planned  for  FY89 .  8 

3.2  Work  Planned  for  FY90 .  8 

4.0  RESULTS  OBTAINED  IN  FY89 .  9 

4.1  Literature  Study .  9 

4.2  PDHOLE  Code  Modification .  9 

4.2.1  Ply  Failure  Criterion .  10 

4.2.2  Stiffness  Degradation  Model .  11 

4.2.3  Programming  Implementation  .  11 

REFERENCES .  13 


ft 


Accession  For  ^  [ 

HI IS  GRAAI 

Sr 

DTIC  TAB 

□ 

Unannounoed 

□ 

Just  If ioatlou- 

b4.  "t-  button/ 


Avt,. lability  Codas 


Olat 


fc! 


lAvall  and/or 
Special 


i 


1.0  INTRODUCTION 


With  the  advances  of  the  composite  materials  technology,  the  use  of  high- 
performance  laminated  composites  for  naval  structures  has  been  substantially 
increased  in  recent  years.  For  examples,  E-glass/epoxy  and  graphite/epoxy  composite 
laminates,  respectively,  have  been  used  by  Naval  Ocean  Systems  Center  to  build 
transducers  and  advanced  unmanned  deep-ocean  search  systems  (AUSS).  In  designing 
these  laminated  composite  structures,  it  is  vitally  important  to  have  a  model  that  can 
accurately  predict  the  strengths  of  the  structures  under  service  conditions.  One  such 
model,  the  progressive  failure  model,  can  take  into  account  the  progressive  damage 
and  fracture  process  of  composite  laminates.  The  model  is  usually  incorporated  into  a 
finite  element  structural  analysis  code  for  composite  structures  due  to  the  complexi¬ 
ties  of  the  structural  problems  encountered. 

An  extensive  literature  review  was  conducted  in  FY88  (via  an  independent 
research  (IR)  program)  on  various  progressive  failure  modeling  schemes  (references  1 
to  7).  Correlations  of  the  analytical  results  predicted  by  these  progressive  failure 
models  with  the  experimental  data  on  composite  specimens  with  complicated  geome¬ 
tries  (e.g.,  hole,  slit,  crack,  etc.)  and  fractured  in  the  in-plane  modes  (e.g.,  fiber  break¬ 
age,  matrix  cracking,  etc.)  were  found  not  totally  satisfactory.  The  ply  failure  crite¬ 
rion  used  in  these  models  may  be  the  prime  reason  for  this.  Meanwhile,  an  improved 
multiaxial  failure  criterion,  the  piecewise  quadratic  strength  tensor  criterion  (refer¬ 
ences  8  and  9),  had  been  recently  developed.  Clearly,  the  introduction  of  the  improved 
ply  failure  criterion  into  the  progressive  failure  model  should  certainly  improve  the 
accuracy  in  predicting  the  strength  of  laminated  composite  structure.  In  light  of  this, 
it  was  proposed  to  use  the  improved  multiaxial  failure  criterion  to  implement  a  pro¬ 
gressive  failure  model  in  a  finite  element  code  for  predicting  the  strength  of  lami¬ 
nated  composite  structure  under  static  loads  and  fractured  in  the  in-plane  failure 
modes. 


The  next  section  provides  the  background  information  for  the  progressive  fail 
ure  model  implementation.  Section  3  discusses  the  implementation  approach 
resulting  from  FY88’s  IR  effort,  formulating  the  work  plan  for  FY89  and  FY90.  Sec¬ 
tion  4  highlights  the  intermediate  results  obtained  in  FY89. 
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2.0  BACKGROUND 


For  background  information,  this  section  briefly  reviews  the  concepts  for  com¬ 
posite  material  stress  analysis  (references  10,  11,  and  12),  the  piecewise  quadratic 
strength  tensor  failure  criterion  (references  8  and  9),  and  progressive  failure  model 
(references  1  to  7).  For  more  detailed  account  of  these  concepts,  including  the  nota¬ 
tion  conventions,  one  may  refer  to  the  references  cited  above. 

2.1  COMPOSITE  STRESS  ANALYSIS 

A  lamina  (ply),  as  shown  in  figure  1(a),  is  a  flat  (sometimes  curved  as  in  a 
shell)  arrangement  of  unidirectional  fibers  or  woven  fibers  in  a  matrix.  A  laminate 
(shown  in  figure  1(b))  is  a  stack  of  laminae  with  various  orientations  of  fiber  direc¬ 
tions  (called  ply  orientations)  in  the  laminae. 

PLY 

ORIENTATION 

ANGLE 


Figure  1.  Laminated  composite. 

Because  of  the  inherent  heterogeneous  nature  of  fiber-reinforced  composite 
materials,  they  are  conveniently  studied  from  two  points  of  view:  micromechanics  and 
macromechanics.  In  micromechanics,  composite  material  behavior  is  studied  by 
examining  the  interaction  of  the  constituent  fibers  and  matrixes  on  a  microscopic 
scale.  In  macromechanics,  composite  material  behavior  is  studied  by  presuming  the 
material  to  be  homogeneous  and  detecting  the  effects  of  the  constituent  materials 
only  as  averaged  apparent  properties  of  the  composite.  Use  of  both  the  concepts  of 
micromechanics  and  macromechanics  allows  the  tailoring  of  a  composite  material  to 
meet  a  particular  structural  requirement  with  little  waste  of  material  capability.  The 
ability  to  tailor  a  composite  material  to  its  job  is  one  of  the  biggest  advantages  that 
composites  have  over  metallic  or  plastic  structures. 

Figure  2  depicts  the  computational  sequences  for  composite  stress  analysis, 
including  those  for  laminated  plate/shell  material  properties  determination  and  ply 
strains/stresses  calculation.  First,  the  micromechanics  properties  (matrix  Youngs 
modulus  and  Poisson's  ratio,  Em  and  vm ;  fiber  Young's  modulus  and  Poisson’s  ratio, 
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Ef  and  Vf ;  and  ply  orientation  angle,  e)  of  each  ply  in  the  laminate  under  considera¬ 
tion  are  homogenized  to  obtain  the  ply  apparent  or  macromechanics  properties  (ply 
(reduced)  stiffness  matrix,  Q ,  whose  components  are  computed  based  on  the  Young's 
modulus  in  the  fiber  direction,  £t ;  the  Young’s  modulus  in  the  direction  transverse  to 
the  fibers,  £2;  the  in-plane  shear  modulus,  Gn;  the  irayor  Poisson’s  ratio,  v12,  etc.). 
Second,  the  macromechanics  properties  of  all  the  plies  are  integrated  through  the 
thickness  of  the  laminate  by  means  of  Kirchhoff  type  of  assumptions  (e.g.,  linear 
strain  variation  through  the  thickness)  to  obtain  the  laminate  (plate/shell)  stiffnesses. 


Homogenization  Kirchhoff 

Process  Type  Assumptions 


(a)  Material  properties  determination 


{ 


2  =  fi  £  • 


(b)  Ply  strains/stresses  calculation 


Figure  2.  Composite  stress  analysis. 


With  the  laminate  material  properties  so  determined  and  boundary  conditions 
(loads  and/or  displacements)  applied,  the  associated  laminated  plate  or  shell  bound¬ 
ary-valued  problem  can  be  solved  by  the  conventional  plate/shell  theories.  The  solu¬ 
tions  contain  reference-surface  (usually  midsurface  of  the  laminated  plate/shell)  defor¬ 
mations:  strains,  £0,  and  curvatures,  £.  Finally,  the  strains,  £,  and  stresses,  £,  for  each 
ply  can  be  computed  by 


Z.-£o  +  z&  . 

(1) 

Otj 

11 

bl 

(2) 

where  2  is  the  coordinate  dong  the  thickness  direction. 
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2.2  THE  PIECEWISE  QUADRATIC  STRENGTH  TENSOR  FAILURE 
CRITERION 


For  a  general  anisotropic  composite,  the  piecewise  quadratic  strength  tensor 
failure  criterion  can  be  written  as1 

/(a*)  =  FiOi  +  FijOiOj  +  HjOj  |  Hjoj  \  =  1.  (i,j,  k  =  1 . 6),  (3) 

where/  is  a  scalar  function,  o*  is  the  contracted  notation  of  the  second-rank  stress 
tensor,2  F,  and  Hi  are  strength  tensors  of  rank  two,  and  Fy  is  a  strength  tensor  of 
rank  four.  In  addition,  constraints  (usually  referred  to  as  the  stability  conditions) 
must  be  imposed  on  the  strength  tensors  Fy  and  Hi  to  ensure  that  the  material 
strength  is  finite  in  all  directions.  More  specifically,  these  constraints  are 

(  Fij  +  HiHj  )oiOj  >  0  (4) 

for  all  stress  points  07  in  the  half  space 

Hi  oi  a  0,  (5) 


and 


(  Fij  -  HiHj  )otOj  >  0 


(6) 


for  all  stress  points  a,  in  the  other  half  space 


Hi  oi  <  0  . 


(7) 


Geometrically,  the  failure  surface  represented  by  equation  (3),  with  the 
strength  tensors  Fy  and  //;  satisfying  the  stability  conditions  given  by  equations  (4) 
and  (6),  is  a  piecewise  ellipsoid  in  the  six-dimensional  stress  space,  which  consists  of 
two  ellipsoids  in  the  two  half  spaces  defined  by  equations  (5)  and  (7).  Hence,  the  fail¬ 
ure  criterion  represented  by  equation  (3)  has  been  referred  to  as  the  piecewise  quad¬ 
ratic  failure  criterion. 

The  above  results  hold  for  a  general  anisotropic  solid.  They  have  been  reduced 
to  an  orthotropic  material,  a  transversely  isotropic  material,  and  an  isotropic  mate¬ 
rial  for  application  to  composites  with  various  material  symmetries.  The  reduced 


'Unless  otherwise  indicated,  the  usual  summation  convention  over  a  repeated  index  is  used  throughout  this  report. 

^With  reference  to  a  rectangular  Cartesian  coordinate  system  (i.e.,  xyi  or,  equivalently,  *1*2X3  system): 

Oi  =  Ox>  Oi  =  Oy,  O3  =  Ot,  O4  =  Oxy,  Oi  =  Oy?,  Ot  -  Ogx  • 
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results,  including  explicit  expressions  for  the  strength  criteria  and  the  restrictions 
imposed  on  the  components  of  the  strength  tensors  occurring  in  these  criteria,  can  be 
found  in  references  8  and  9. 

2.3  PROGRESSIVE  FAILURE  MODEL 

Estimates  for  ultimate  strengths  of  fiber-reinforced  laminates  based  on  first 
ply  failure  concept  are  often  highly  conservative  because  first  ply  failure  usually  does 
not  result  in  the  total  failure  of  the  laminate.  To  make  reasonable  strength  prediction 
at  ultimate  failure  of  the  laminate,  account  has  to  be  taken  into  the  progressive  dam¬ 
age  of  the  laminate  and  the  resulting  stress  redistribution  that  can  occur  due  to  such 
damage. 

To  show  how  these  can  be  taken  into  account  by  a  progressive  failure  model, 
the  numerical  procedures  for  the  model  incorporated  in  a  finite  element  code  are 
given  in  figure  3.  These  procedures  are  performed  according  to  the  following  steps: 

1.  Increase  the  applied  load/displacement  by  a  small  increment  and  perform 
composite  stress  analysis  as  described  in  subsection  2.1  to  obtain  strains/ 
stresses  for  each  element  in  each  ply  of  the  laminate. 

2.  Access  damage  by  comparing  the  computed  stresses/strains  with  the  ply 
(and/or  interply)  failure  criterion.  Return  to  the  first  step  if  no  damage  is 
found.  Continue  to  the  following  steps  if  damage  occurs. 

3.  Degrade  (modify)  ply  or  interply  macromechanics  properties  (e.g.,  ply  stiff¬ 
ness  Q)  for  the  damaged  (failed)  elements. 

4.  Perform  composite  stress  analysis  for  stress  redistribution. 

5.  Repeat  steps  2  to  4  until  no  more  damages  can  be  found  at  the  given  load/ 
displacement  increment. 

6.  Repeat  steps  1  to  5  until  equilibrium  can  no  longer  be  achieved  or  no  more 
load  can  be  sustained. 

From  figure  3,  it  should  be  clear  that  the  following  modules  are  essential  for  a 
progressive  failure  model  implemented  in  a  finite  element  code: 

a.  Progressive  failure  modeling  driver  for  load  increment  loop  and  stress 
redistribution  loop  at  a  given  load  increment, 

b.  Ply  failure  criterion, 

c.  Stiffness  degradation  model, 

d.  Damage  propagation  output. 

Before  leaving  this  section,  note  that  various  schemes  of  progressive  failure 
modeling  differed  from  one  another  primarily  in  the  ply  failure  criterion  as  well  as 
the  degradation  model  of  the  stiffness  of  the  laminate  under  failure  of  one  or  more 
plies. 
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Figure  3.  Progressive  failure  model. 


3.0  DEVELOPMENT  PLAN 


To  implement  the  proposed  progressive  failure  model  into  a  general-purpose 
finite  element  code,  several  such  codes  for  laminate  composite  structural  analysis 
were  considered  in  FY88.  These  included: 

a.  CODSTRAN  —  Composite  Durability  STRuctural  ANalysis  (reference  13), 

b.  ABAQUS  (reference  14), 

c.  MSC/NASTRAN  —  MacNeal-Schwendler  Corporation/NAsa  STRuctural 
ANalysis  (reference  15), 

d.  COSMIC/NASTRAN  —  COmputer  Software  Management  and  Information 
Center/NAsa  STRuctural  ANalysis  (reference  16). 

Among  these  codes,  COSMIC/NASTRAN  was  found  to  be  the  only  code  that 
encompasses  all  of  the  following  features: 

a.  Availability  of  the  source  codes  for  programming  implementation, 

b.  Maturity  of  the  code  developments, 

c.  Frequent  use  in  our  day-to-day  work  at  Naval  Ocean  System  Center. 

Thus,  COSMIC/NASTRAN  was  chosen  for  incorporating  the  proposed  progressive 
failure  model.  This  code  has  been  leased  on  an  annual  basis  and  is  currently  opera¬ 
tional  at  this  center. 

In  addition  to  acquiring  COSMIC/NASTRAN,  the  special-purpose  finite  ele¬ 
ment  program,  PDHOLE,  developed  by  Stanford  University  (reference  7),  has  been 
purchased.  This  program  was  written  specially  for  predicting  the  tensile  strength  of 
a  composite  laminate  strip  containing  an  open  hole.  Due  to  this  special-purpose 
nature,  this  program  is  much  simpler  (shorter)  than  COSMIC/NASTRAN  to  follow. 
Furthermore,  PDHOLE  is  already  equipped  with  the  progressive  failure  modeling 
scheme,  except  with  different  ply  failure  criterion  and  stiffness  degradation  model 
from  what  we  proposed.  In  light  of  this,  the  improved  multiaxial  failure  criterion  as 
well  as  the  associated  stiffness  degradation  model1  can  be  implemented  into  the 
PDHOLE  code  to  gain  failure  model  implementation  experiences.  Moreover,  the  im¬ 
plemented  version  of  the  PDHOLE  code  can  be  used  for  the  subsequent  program  com¬ 
parison  and  check-out  purposes  with  the  implemented  version  of  COSMIC/ 
NASTRAN. 

Thus,  the  required  works  for  implementing  the  two  codes  were  planned  for 
FY89  and  FY90  in  the  following  two  subsections. 


’The  associated  stiffness  degradation  model  will  be  presented  in  subsection  4.2. 
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3.1  WORK  PLANNED  FOR  FY89 


Work  planned  for  FY89  consists  of 

a.  Code-Supporting  Literature  Studies  —  Study  in  detail  all1  the  literature 
directly  relevarJ  cO  both  the  understanding  of  the  two  computer  codes  and 
the  programming  of  the  proposed  progressive  failure  model  into  the  two 
codes. 

b.  PDHOLE  Modifications  —  Understand  the  program  listings  and  program 

the  improved  multiaxial  failure  criterion  and  the  associated  stiffness  deg¬ 
radation  model  into  the  computer  program.  » 

3.2  WORK  PLANNED  FOR  FY90 


Work  planned  for  FY90  consists  of 

a.  COSMIC/NASTRAN  Modifications  —  Run  COSMIC/NASTRAN  with  test 
problems;  study  the  program  listings  for  implementation  purposes; 
implement  a  progressive  failure  modeling  driver  into  the  program;  write 
subroutines  for  the  improved  failure  criterion  and  the  proposed  stiffness 
degradation  model;  and  program  damage  propagation  output  subroutines. 

b.  Code  Verifications  —  Numerically  verify  the  modified  COSMIC/ 
NASTRAN,  using  the  modified  PDHOLE. 


'Except  COSMIC/NASTRAN  program  manuals. 
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4.0  RESULTS  OBTAINED  IN  FY89 


The  work  scheduled  for  FY89  as  described  in  subsection  3.1  has  been  all  com¬ 
pleted.  The  results  are  highlighted  below  in  subsections  4.1  and  4.2. 

4.1  LITERATURE  STUDY 

To  understand  the  PDHOLE  code,  a  series  of  papers  (references  6,  7,  and  17  to 
25)  contributing  to  the  code  development  during  1982  and  1987  have  been  studied  in 
♦  great  depth.  They  provide  insights  into  how  the  following  basic  concepts  have  been 

implemented  into  the  code: 

a.  Progressive  failure  modeling  driver, 

b.  Ply  failure  criterion, 

c.  Stiffness  degradation  model, 

d.  Nonlinear  shear  stress-strain  law  of  lamina, 

e.  Stress  redistribution,  etc. 

Due  to  the  lack  of  good  supporting  documents  of  COSMIC/NASTRAN  and  the 
close  programming  relationship  between  COSMIC/NASTRAN  and  MSC/NASTRAN,1 
understanding  of  COSMIC/NASTRAN  began  by  studying  the  book  MSCINASTRAN 
Primer  (reference  26).  This  book  provided  the  reader  not  only  with  a  description  of 
the  technological  content  of  the  code,  but  also  with  a  description  of  the  NASTRAN 
vocabulary  and  capability  for  structural  analysis.  In  addition  to  the  study  of  this 
book,  a  short  course,  entitled  “DMAP  and  Database  Application  in  MSC/NASTRAN 
Version  66,”  was  taken.  This  course  showed  the  student  how  to  perform  tasks  in 
NASTRAN.  Such  tasks  included  the  creation,  storage,  and  maintenance  of  user- 
written  solution  sequences,  writing  structured  DMAP,  and  operation  of  complex  data¬ 
base  using  the  File  Management  System. 

4.2  PDHOLE  CODE  MODIFICATION 


This  subsection  records  the  results  for  in-plane  ply  failure  criterion  and  the 
associated  stiffness  degradation  model  that  have  been  implemented  in  the  PDHOLE 
code.  The  detailed  derivations  of  these  results  from  the  general  three-dimensional 
theories  are  omitted  here,  but  they  shall  be  given  in  a  forthcoming  final  report  on  the 
progressive  failure  model  development.  Moreover,  this  subsection  also  summarizes 
the  programming  implementation  of  the  original  version  of  the  PDHOLE  code. 


'As  is  well  known,  COSMIC/NASTRAN  and  MSC/NASTRAN  have  common  origins  in  Level  15.5  NASTRAN  (refer¬ 
ence  26). 
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4.2.1  Ply  Failure  Criterion 


In  PDHOLE,  consideration  is  restricted  to  the  plane-stress  state:1 


Oi  *  0 ,02*  0.  On  *  o,  03  =  (Tjj  =  023  =  0, 


(8) 


or 

&=  \  o2 

l  012 

Also,  the  composite  lamina  under  consideration  has  a  material  symmetry  no  more 
general  than  orthotropic.  For  this  material  in  the  plane-stress  state,  the  stress-strain 
relations  are  given  by 


o  =  Q  e. . 


(10) 


where  f  is  the  (ply)  strain  vector  having  three  components: 


and  Q  is  the  (ply)  reduced  stiffness  matrix  given  by 

[Qn  Qn  0  "I 

O12  Qn  0  I  (12) 

0  0  Qu  J  . 

For  an  orthotropic  composite  lamina  in  the  plane-stress  state,  the  general 
multiaxial  failure  criterion  represented  by  equation  (3)  reduces  to2 


f  -F\  0\  +  F2O2  +  ^11  0?  +  2F12  01  02  +  F 22  of  +  Fu  012 

(13) 

+  (  Hi  0\  +  M2  02  )  I  M\  0\  +  H2  02  I  =  1  , 

and  there  are  eight  independent  strength  constants  (Ft ,  F2  \  Fu,  Fn,  F12,  Fu\  Mh  Hx) 
appearing  in  the  ply  failure  criterion.  The  number  of  independent  strength  constants 


'Throughout  the  remainder  of  this  report,  the  1-2  plane  ia  taken  to  be  tangent  to  the  lamina  surface  and  the  3-axis 
is  taken  to  be  along  the  thickness  direction  of  the  lamina  (see  figure  1(a)). 

*In  view  of  the  footnote  on  page  4,  it  can  be  seen  that  0i2  =  0«  . 
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can  be  reduced  to  four  for  either  a  transversely  isotropic  or  an  isotropic  composite 
lamina  in  the  plane-stress  state.  For  either  material,  the  independent  strength  con¬ 
stants  are:  Fx ,  Fn,  Fn;  hu  and  the  following  relations  hold: 

Fi  =  F\  ,  F22  =  Fu  ,  F44  =  2  (  Fu  -  F\2  ),  Hi  -  H 1  .  (14) 

4.2.2  Stiffness  Degradation  Model 

Once  failure  occurs,  the  material  may  undergo  some  degree  of  property  loss  in 
the  damaged  area.  At  the  failed  material  point,  this  is  realized  by  reducing  the  longi¬ 
tudinal  modulus,  £1,  the  tranverse  modulus,  Elt  the  in-plane  shear  modulus,  C12,  and 
the  mqjor  Poisson’s  ratio,  V12,  to  near  zero.  Consequently,  at  the  failure  point,  the 
components  of  the  ply  reduced  stiffness  matrix  Q  can  be  set  as 

Qu  =  Qn  -  Q22  =  Qu  -  10"1  F.\  ,  (15) 

and  all  in-plane  stress  components  are  vanishing: 

0\  —  0\2  —  o2  =  0  .  (16) 


4.2.3  Programming  Implementation 

For  programming  the  proposed  ply  failure  criterion  and  stiffness  degradation 
model  into  the  PDHOLE  code,  a  COMMON  block  named  MP1  has  been  added  to  the 
appropriate  subroutines  in  the  code  and  various  subroutines  have  been  modified. 

The  MP1  COMMON  block  is  given  by: 

COMMON !MP\/ IF AIL  ,  F\  ,  F2  ,  F\\  ,  F\2  ,  F22  ,  F44  ,  HI  ,  H2  ,  (17) 


where 


IFAIL  =  Failure  Criterion  Option 

=  0:  Use  the  failure  criterion  originally  programmed  in 

the  code, 

1 :  Use  the  proposed  piecewise  quadratic  strength 
tensor  failure  criterion. 

FI  ,  F2  ,  FI  1  ,  F12  ,  F22  ,  F44  ,  HI  ,  H2 

=  FORTRAN  variable  names  for  the  strength  constants 
input:  Fi  ,  Fj  ,  Fn  ,  F12  ,  F22  1  F44  ,  H\  ,  Hi . 


Among  those  which  have  been  extensively  modified  are  the  following  subrou¬ 
tines: 


MAIN 

Main  program. 

INPUT 

Input  material  properties, 

CFAIL 

Ply  failure  criterion, 

PMATRL  - 

Stiffness  degradation  model, 

STRESS  - 

Stress  calculation, 

RDSTR 

Stress  redistribution. 
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